Sublattice synchronization of chaotic networks with delayed couplings 
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Synchronization of chaotic units coupled by their time delayed variables are investigated analyti- 
cally. A new type of cooperative behavior is found: sublattice synchronization. Although the units 
of one sublattice are not directly coupled to each other, they completely synchronize without time 
delay. The chaotic trajectories of different sublattices are only weakly correlated but not related by 
generalized synchronization. Nevertheless, the trajectory of one sublattice is predictable from the 
complete trajectory of the other one. The spectra of Lyapunov exponents are calculated analytically 
in the limit of infinite delay times, and phase diagrams are derived for different topologies. 
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Two identical chaotic systems which are coupled to 
each other can synchronize to a common chaotic trajec- 
tory. Both of the systems are chaotic, but after syn- 
chronization the two chaotic trajectories are locked to 
each other [l|, This phenomenon has attracted a 
lot of research, partly because chaos synchronization has 
the potential to be applied for novel secure communi- 
cation systems 0, Hjj- I n fact, synchronization and bit 
exchange with chaotic semiconductor lasers has recently 
been demonstrated over a distance of 120 km in a public 
fiber-optic network [||. 

Typically, the coupling between chaotic units has a 
time delay due to transmission of the exchanged signal. 
Nevertheless, two chaotic systems can synchronize with- 
out delay, isochronically, although the delay time may 
be very long. This counterintuitive phenomenon has re- 
cently been demonstrated with chaotic semiconductor 
lasers @, 0], and it is discussed in the context of cor- 
responding observations on correlated neuronal activity 

In this Letter we give an analytic solution of small 
networks of chaotic units coupled by their time-delayed 
variables. Motivated by chaotic lasers, we consider the 
competition between delayed self-feedback and coupling. 
For regular networks which can be decomposed into two 
identical sublattices - e.g. a ring with N = 4 units in 
Fig- IT|b) ~ a new kind of synchronization is observed: 
The units (A, C) as well as the units (B, D) are com- 
pletely synchronized. Even for arbitrarily long delay 
times, the synchronization is isochronical, i.e. the trajec- 
tories of units A and C are identical without any time- 
shift. Although synchronization of (A, C) is caused by 
the trajectory of (B, D) and vice versa, the trajectories 
of A and B are not synchronized. Actually, we do not 
even find generalized synchronization between A and B. 
The cross-correlations are symmetric, there is neither a 
leader nor a laggard for the two sublattices. Each unit 
is passive, i.e. it has negative Lyapunov exponents when 
it is isolated and driven by an external signal. Neverthe- 
less, the mutual interaction of passive elements leads to 



high-dimensional chaos with synchronized units in each 
sublattice. 
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FIG. 1: Small networks of chaotic units. Double lines signify 
bidirectional couplings whereas arrows show unidirectional 
couplings. It turns out that each ring with an even num- 
ber N of units has Lyapunov exponents which are identical 
to a chain with N/2 + 1 units. 



We have found sublattice synchronization in several 
chaotic networks, in particular for rate equations describ- 
ing chaotic semiconductor lasers 
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In this Letter we 
analyze sublattice synchronization using a simpler sys- 
tem, namely chaotic Bernoulli maps with delay, which 
allows analytical as well as comprehensive numerical in- 
vestigations. In the limit of long delay times, the spec- 
tra of Lyapunov exponents, which are calculated analyti- 
cally, determine the phase boundaries between sublattice 
synchronization, complete synchronization and complete 
chaos. 

Two mechanisms are included in our model: 1. De- 
layed feedback and 2. delayed mutual coupling. The 
first mechanism can generate high dimensional chaos for 
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each single unit whereas the second one synchronizes the 
network. The phenomena caused by these two mecha- 
nisms are captured by a network of iterated maps, given 
by the following equations, j, k = 1, . . . , N: 



x J t 



(l- £ )/(4_ 1 ) + 



E 

k 



F Jtk f(xt T ) (1) 



Each node j of the network contains a variable x J t £ 
[0, 1] which is iterated in time t. Fj t k is the weighted 
adjacency matrix which has a component 1/zj if node j 
is driven by node k and otherwise, and Zj is the number 
of connections to node j. 

Analytic calculations are possible for the Bernoulli 
shift map f(x) = (ax) mod 1 which is chaotic for a > 1 
r is the delay time, which, for simplicity, is identical 
for the feedback as well as for the coupling. The param- 
eter e determines the strength of the delay terms while 
k determines the relative strength of the self-feedback 
compared to the mutual coupling terms. 

Now let us discuss stationary solutions of IJJ). Com- 



plete synchronization, 



i 



x\ for all (j, k), is a solu- 



tion of |JJ . In addition, several kinds of incomplete syn- 
chronization are solutions, as well. For example, sublat- 
tice synchronization, x\ — x\ for all nodes (j, k) in each 
sublattice, is a solution, too, if the network can be de- 
composed into two interconnecting sublattices, as for the 
chains and rings in Fig. [TJ But are these solutions sta- 
ble? Their stability has to be determined from equation 
([TJ). Let 8x 3 t be a small deviation from any given trajec- 
tory x\ . The linearized equations ([TJ are solved with the 
ansatz Sx J t — c Vj, where the value A = In \c\ is the corre- 
sponding Lyapunov exponents. We find that the vectors 
(Vj)j=i,..,N are eigenvectors of the adjacency matrix F, 
and their corresponding eigenvalues fi q ,q = 0, ...,N — 1 
determine the Lyapunov exponents by 



(1 — e)ac T 1 + sua 
e(l — k)q 



(2) 



For each eigenvector, indexed by q — 0, N — 1, this 
polynomial equation |(5J) for c has r complex solutions 
Cq,m , tci — 1 , . . . , t giving a total of Nt Lyapunov expo- 
nents which are plotted in Fig. [5] The eigenvectors of 
F determine the N directions of perturbation v q j which 
shrink or explode with the corresponding values of c. 

Here we are interested in long delay times, r — > oo. It 
turns out that for each eigenvalue fj, q , q — 0, . . . , N — 1, 
there is one Lyapunov exponent of order one, and r — 1 
exponents of order 1/r (see Fig. [2]). The large exponent 
is given by = ln[(l — e)a] and merges with the band 
of the other r — 1 exponents at the critical point 
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FIG. 2: Lyapunov spectra for the different modes q for TV = 
4, a = 1.5, t = 40, k = 1/4. At e = 4/9 the (q = 1) mode 
becomes stable (— > sublattice synchronization). At e = 2/3 
the (q = 2) mode, (+1, — 1, +1, — 1), becomes stable, too (— » 
complete synchronization) . See also Fig. [3] 



Consequently, for e < e c , stable synchronization does not 
exist. 
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FIG. 3: Phase diagram for a = 1.5. See text for assignment 
of the different regions to complete and sublattice synchro- 
nization. 



The eigenvalues of the adjacency matrix F determine 
the stability of synchronized trajectories according to (JSJ, 
and they depend on the topology of the lattice. Let us 
consider the rings and chains depicted in Fig. [TJ We find 
that the phase diagram of all six topologies is given by 
Fig. El To see this, let us consider a ring of N mutu- 
ally coupled chaotic units with self-feedback, described 
by UJ. Figure Q] shows N — 2, 4 and 6. By symmetry, 
the eigenvectors v q .j and eigenvalues fi q of the adjacency 
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matrix F are the Fourier modes, (q = 0, N — 1) 



l Q:> 



exp(2mj q/N), \x q = cos J . (4) 



These eigenvalues determine the Lyapunov spectrum by 
©. Using the mehtods of Ref. [l2j . for large values of r 
the maximal Lyapunov exponent is calculated as 



A" lax = - In 

q T 



(It 
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1 -a(l -e) 



(e>e c ). (5) 



The mode (q = 0) with the corresponding eigenvec- 
tor (1, 1,1,.. .) is always unstable, AJ^fg > ^, hence the 
system is chaotic for all parameters e and k. In fact, 
the system is in a state of high-dimensional chaos (hy- 
perchaos), since the Kaplan- Yorke dimension increases 
linearly with the delay time r. 

The superpositions of the modes (q = l,q = N — 1), 
which have degenerate eigenvalues, are perturbations 
which destroy any structure of synchronization. They 
are unstable for 



k > K, 



1 



9=1 



as 



[1 



COS I 



N 



)] 



(6) 



An analysis of shows that for even values of N, 
the only perturbation which is relevant in the region k < 
Kg— i is the perturbation of the mode (q = N/2). The 
network is unstable against this mode for 



K < K q= N/2 — 



a - 1 

2as 



(7) 



Since this mode corresponds to the eigenvector 
(+1, — 1, +1, — 1, . . . , — 1), it can destroy only complete 
but not sublattice synchronization. 

Equations ([6]) and (O immediately determine the 
phase diagram of the rings with N units in the (k, e) 
plane, Fig. O The ring with N = 4 units, Fig. QJb), 
is completely synchronized in region II (= Ha plus lib). 
But complete synchronization is cannot exist in region 
III, because there the q = N/2 mode is unstable. In this 
region the system is synchronized in a sublattice con- 
figuration: the trajectory of A is completely identical, 
without delay, with the one of C, and correspondingly, B 
is identical to D. 

When increasing the number N of units the (q = 1) 
line of Fig. [3] shrinks to the lower right corner. For N = 6 
wc find complete synchronization in the region lib, while 
sublattice synchronization occurs in region Illb. In the 
latter region the trajectories of (A, C, F) as well as (B, 
D, G) are identical, but the two different trajectories do 
not synchronize. 

Finally, for large system size N synchronization dis- 
appears completely. Complete synchronization exists up 
to 



N < 



2tt 



arccos 



(8) 



while sublattice synchronization remains for larger sizes 
up to 



N < 



2n 



(9) 



Consequently, even for arbitrary large rings one obtains 
regions of complete and sublattice synchronization, if the 
chaos is sufficiently weak. Asymptotically, sublattice syn- 
chronization is stable if a is smaller than a ~ l+2(7r/jV) 2 . 

In our model the rings with an even number N of 
chaotic units are equivalent to chains with N/2 + 1 units, 
i.e. the corresponding rings have identical, but degener- 
ate, Lyapunov exponents. Hence the phase diagram of 
Fig. [3] describes synchronization of the chains Fig. [lja), 
(d) and (e), as well. Regions II and III also describe the 
infinite chain with directed bonds of Fig. QJf). If the left 
unit A is just a mirror , corresponding to a self-feedback 
of B with k = 1, the chain is in state of complete syn- 
chronization in these regions. If, however, A is a chaotic 
unit, the infinite chain switches to a state of complete 
synchronization in region II and sublattice synchroniza- 
tion in region III. 

These findings are not specific to rings and chains, 
since we found sublattice synchronization in square lat- 
tices with periodic boundaries (with even side lengths) 
as well as for free boundaries (with even or odd side 
length). In addition, we found sublattice synchroniza- 
tion for a small triangular lattices with three sublattices. 
Each sublattice synchronizes completely, without time 
shifts, whereas the three trajectories have only weak cor- 
relations. 

Synchronization of the units of one sublattice is relayed 
by the chaotic trajectory of the other one. Hence one 
would expect some relation between these two chaotic 
trajectories (or three for the triangular lattice). 

The cross correlations between the two sublattices are 
shown in Fig. [U For sublattice synchronization the cen- 
tral peak, a reminder of complete synchronization, still 
exists, in addition to correlations shifted by multiples of 
the delay time r. Only when self-feedback is switched 
off, k = 0, the isochronal correlations disappear. The 
correlations are completely symmetric in time shift ±A, 
there is no symmetry breaking to a leader/laggard sce- 
nario. In addition, when crossing the phase boundary be- 
tween sublattice synchronization, complete synchroniza- 
tion and complete chaos, the central peak of the cross- 
correlation changes discontinuously, similar to a subcrit- 
ical bifurcation. For large, but finite observation times, 
our simulations even show hysteresis effects. 

Apart from these weak correlations, is there any func- 
tional relation between the two different trajectories? In 
particular, are the trajectories in a state of generalized 
synchronization [l3|]? For our model, generalized syn- 
chronization means that there exists a function between 
the sequence (xf, . . . , xf_ T ) of A and [xf , . . . , xf_ T ) of 
B. This function may be fractal [l[ or may have several 
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FIG. 4: Cross-correlation C between At and B t +A for a = 
1.5, t = 40, £ = 0.7. (a): Complete Synchronization, (b,c) 
sublattice synchronization. 



will relax to a unique trajectory. Since we find sublattice 
synchronization only for passive units, we can predict the 
trajectory of B from the knowledge of the complete tra- 
jectory of A after some transient time. Note that for two 
units we find complete synchronization for active units, 
as well, but in this case, region I, prediction is trivial. 

Sublattice synchronization has not yet been found in 
experiments. However, a chain of three lasers Fig. QJd), 
has been investigated with semiconductor lasers without 
feedback, k — :7j]. Synchronization of the outer lasers 
has been observed, being relayed by a different chaotic 
trajectory of the internal laser, in agreement with our 
model. But there are differences to our results: In our 
model, the cross correlation is symmetric in time shift, 
there is neither a spontaneous symmetry breaking nor a 
leader/laggard scenario [15J, and there is no generalized 
synchronization between the synchronized and the con- 
necting units 

We would like to thank Jorg Reichardt, Andreas Rut- 
tor and Evi Kopelowitz for useful discussions. 



branches, which usually makes an analysis difficult. But 
at least for t = 1 of r = 2 the relation between these two 
vectors can be easily visualized. Fig. \5\ shows that the 
trajectories are not related by generalized synchroniza- 
tion in the phase of sublattice synchronization. 
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FIG. 5: Two trajectories of the sublattice synchronization 
phase, for the networks of Fig. [T]with r = 1, e = 0.8, k — 
0, a — 3.0. Obviously, there is no functional relation between 
the two trajectories. 

However, there is a related question: Can the trajec- 
tory of sublattice B be predicted from the knowledge of 
the complete trajectory of A 14|? This question has a 
positive answer for the configurations of Fig. [TJ In re- 
gions II and III all units are passive. This means, if we 
omit the drive, the last term of ([I]), each unit has nega- 
tive Lyapunov exponents. Hence if this unit is driven, it 
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